First principles study of ceramic materials (IVB group carbides) under ultrafast laser irradiation
He Nan-Lin1, Cheng Xin-Lu1, 2, †, Zhang Hong1, 2, 3, Yan Gai-Qin3, Zhang Jia1
Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
Key Laboratory of High Energy Density Physics and Technology of Ministry of Education, Sichuan University, Chengdu 610065, China
College of Physical Science and Technology, Sichuan University, Chengdu 610065, China

 

† Corresponding author. E-mail: chengxl@scu.edu.cn

Abstract

Group IVB carbides have been applied in extreme aerospace environments as hard ceramic coatings; ZrC is being considered as a replacement for SiC in nuclear reactors. Therefore, a thorough understanding of the laser irradiation response of group IVB carbides is of clear significance. However, the existing knowledge on the fundamental properties of IVB group carbides is limited and insufficient with regard to both irradiated and non-irradiated characteristics. We investigate the effect of ultrafast laser irradiation on the lattice stability of ceramic materials (IVB group carbides) using the density functional perturbation theory (DFPT). The calculated phonon frequencies of TiC and ZrC at the ground state are in good agreement with previous calculations and experimental values. The phonon frequencies of IVB group carbides are positive, even though the electronic temperature reached 5 eV. Thus, IVB group carbides are more stable under ultrafast laser irradiation, which has greater benefits in nuclear and aeronautical applications compared to metals (W, Na), semimetals (Bi), and semiconductors (Si, SiC). The thermodynamic properties of ZrC are calculated as functions of their lattice temperature at different electronic temperatures. The elastic shear constants of IVB group carbides satisfy the Born stability criteria at Te = 5 eV. In addition, a comparison of the predicted melting temperatures of IVB group carbides, reveal that HfC is better suited for extreme high-temperature environments.

1. Introduction

Ceramic materials (IVB group carbides) constitute different degrees of the three characteristic bonding forms: ionic, covalent, and metallic. These materials show good corrosion/oxidation resistance, very high melting point, extreme hardness, and other covalent properties. In addition, they exhibit strikingly high electrical/thermal conductivities compared with those of pure transition metals. Moreover, they have a NaCl-type structure that is associated with ionic bonding.[15] By exploiting their favorable physical and chemical properties, on the one hand, the transition metal carbide refractory ceramics can be candidate materials for use in high-temperature nuclear applications. In particular, ZrC is being considered as a coating material for the tri-isotropic coated nuclear fuel used in high temperature reactors, replacing or in addition to the currently used SiC.[68] Shen et al.[9] studied the lattice stability of β-SiC under ultrafast intense laser irradiation. To the best of our knowledge, there are no published reports on ZrC. On the other hand, transition metal carbides have numerous applications in the aeronautical industry.[10] Amir et al.[11] reported that the utilization of ceramic engine for the future aerospace applications may result in enhanced performance characteristics and reduced operational costs. Structural materials are typically exposed to various types of irradiation, including ultrafast particle and laser irradiation, in both nuclear applications and the aeronautical industry. The effects of fast neutron irradiation (1−10×1025 N/m2 and 635–1480 °C) on the properties of high-purity zone refined zirconium carbide were measured by Snead et al.,[12] including the lattice parameter, hardness, elastic modulus, thermal conductivity, and indentation fracture toughness. They discovered that this ceramic is quite stable under fast neutron irradiation in the temperature and dose range. The processes of damage evolution in TiC crystals irradiated with hydrogen and deuterium ions at low temperatures were examined in situ with an electron microscope equipped with an ionic accelerator. Amorphization was confirmed upon hydrogen ion irradiation, while no amorphization occurred during deuterium ion irradiation.[13] The results of a previous study on the effect of ion irradiation on the microstructure stability of GFR ceramics (ZrC, ZrN, TiC, TiN, and SiC irradiated with 1 MeV Kr-Ions to 10 and 70 dpa at 800 °C) was reported.[14] Huang et al.[15] published results on the damage evolution of ZrCX (where x ranges from 0.9 to 1.2) under proton irradiation at 800 °C, and discussed the irradiation-induced defects, such as density of dislocation loops, at different stoichiometries and doses. Changes in the physical properties under a range of irradiation and high temperature conditions are very important in nuclear and aeronautical applications. However, previously reported results on laser irradiation research on transition metal carbide refractory ceramics are limited to microsecond laser pulse interactions. There are few published calculations and experiments which specially focus on the effect of femtoseconds ultrafast intense laser irradiation on IVB group carbides.

Irradiating a target material using an ultrafast laser pulse (∼100 fs) can excite a large number of electrons from the valence band to the conduction band. The electron-electron collisions result in the electronic energies approaching a Fermi-Dirac distribution with a well-defined electronic temperature Te of approximately 104 K in a few tens of femtoseconds. However, the ion subsystem remains close to its initial temperature. Such a high electronic temperature would dramatically change the electronic charge distribution and the inter-atomic potential, inducing different physical effects on the target materials, such as phonon squeezing,[16] solid-solid transitions,[1720] and non-thermal melting.[2123] Herein, we perform calculations for the lattice dynamic, thermodynamic and elastic properties, and the predicted melting temperatures of IVB group carbides at different electronic temperatures.

The present investigation is organized as follows. In section 2, the computational method and technical details are described. Section 3 discusses the lattice dynamic, thermodynamic and elastic properties, and the predicted melting temperatures of IVB group carbides at different electronic temperatures. Finally, conclusions are summarized in Section 4.

2. Computational method and technical details

The calculations were performed with the Vienna ab initio simulation package (VASP) code[2426] based on density functional theory (DFT). The projector-augmented wave (PAW)[27] method was employed to describe the electron-ion interactions. For all the geometrical optimizations and the calculations of properties, we adopted a face-centered cubic (fcc, space group , 225) structure. The exchange–correlation effects are mainly taken into account using the Perdew–Burke–Ernzerhof (PBE) function[28] of the generalized gradient approximation (GGA).[29] For TiC, ZrC, and HfC, the valence electrons configurations are 3d24s2 for Ti, 4d25s2 for Zr, 5d26s2 for Hf, and 2s22p2 for C, while the more tightly bound electrons are represented as core electrons. The plane-wave basis sets with the energy cutoff are set to 700 eV for TiC and 500 eV for ZrC and HfC. An 18×18×18 Monkhorst–Pack mesh of k points in the first Brillouin zone (BZ) was sampled for geometrical optimization at the different electronic temperatures. All parameters were tested for convergence. Moreover, suitable bands were employed to ensure the availability of enough states for electrons to occupy even at high electronic temperatures.

Phonon frequencies were obtained by inter-atomic force constants in the real space[30,31] within density functional perturbation theory (DFPT)[32] at different electronic temperatures. In the calculation, we use the 2×2×2 supercell and the 6×6×6k mesh to acquire the real space force constants. In addition, the thermodynamic functions (such as the phonon free energy F, the phonon internal energy E, the phonon entropy S, and the phonon heat capacity CV) can be calculated from the phonon frequencies w = (q, l) and phonon density of states g(w) through the following functions.[3234]

where N is the number of primitive unit cells, n the number of atoms per primitive unit, h the Planck constant, kB the Boltzmann constant.

At the same time, the total elastic moduli matrix was studied using the stress-strain method based on DFT with the 4×4×4k mesh. For cubic crystals, the Voigt average shear modulus was used to analyze the elastic stability of the structure, where C44 and are the two shear moduli.[35] From the results of the elastic moduli, the Debye temperature ΘD at different electronic temperatures can be estimated using the following equation[36]

where vm is the velocity of sound, k the Boltzmann constant, ρ the density, M the molecular weight, NA the Avogadro number, and n the number of atoms in the molecule. The relationship between the Debye temperature ΘD and the melting temperature Tm is assumed by the Lindermann melting criterion:[37] , where A depends on the density and the atom mass. We can therefore predict the melting temperature at different electronic temperatures utilizing the following relation: .

3. Results and discussion
3.1. Lattice optimization at different electronic temperatures

In order to investigate the effect of ultrafast laser irradiation on the lattice parameters, we first calculated the lattice parameters for the ground state. The values of TiC, ZrC and HfC are 4.3375 Å, 4.7239 Å, and 4.647 Å, respectively, which are very close to the experimental values[38] of 4.332 Å, 4.692 Å, and 4.639 Å, respectively. A comparison with experimental data revealed an overestimation of the equilibrium lattice parameters by 0.01% for TiC, 0.06% for ZrC and 0.02% for HfC. These results are in agreement by GGA standards. We then added the different electronic temperatures for the systems and obtain the optimized lattice parameters of the IVB group carbides. Figure 1 shows that the lattice parameters of the IVB group carbides clearly increased with an elevation of the electronic temperature. In addition, both TiC and ZrC exhibit approximately the same trend. As is well known, the systems are in mechanical equilibrium in the ground state. However, with the elevation of the electronic temperature, the previous equilibrium of the systems is destroyed by the electronic excitation. The inter-atomic attractive interactive would be weakened while the repulsive interaction is almost not affected by the excited electrons. To achieve a new equilibrium, the lattice volume expands, namely, the lattice parameters increases with an increase in the electronic temperature.

Fig. 1. (color online) The lattice parameters of IVB group carbides at different electronic temperatures.
3.2. Lattice dynamical properties at different electronic temperatures

We calculated the phonon dispersion curves of the IVB group carbides in fcc as a function of electronic temperature to investigate the lattice dynamic stability under ultrafast laser irradiation. The different imaginary phonon frequencies at some points of the Brillouin zone provided different information about lattice dynamic instability.[39] An imaginary frequency of the whole transverse acoustic branch at a high electronic temperature indicates a non-thermal melting transition.[40] An imaginary frequency of the whole longitudinal optical branch at the point of the Brillouin zone indicates coherent optical phonon generation.[41] An imaginary frequency of a branch in the whole Brillouin zone is attributed to disordering of the lattice.[39]

In Table 1, our calculated values of phonon frequencies at several high-symmetry points ( , X, and L) of the Brillouin zone of TiC and ZrC at ground state are shown and compared with previous experimental[42,43] and theoretical values.[44] From Table 1, we can see that our calculated values are in excellent agreement with previous published experimental and theoretical data. This suggest that our results are reliable for the further studies at high electronic temperatures.

Table 1.

The calculated phonon frequencies (cm−1) of TiC and ZrC at ground state and their comparisons with experimental and other theoretical values.

.

The phonon dispersion curves at different electronic temperatures are represented in Fig. 2. The optical and acoustic branches of the phonon spectrum of IVB group carbides holistically decreased with the elevation of the electronic temperature. Partial phonon frequencies begin to include imaginary frequencies, which signify the instability of the lattice at Te = 5.4 eV for TiC, Te = 5.17 eV for ZrC and Te = 5.5 eV for HfC. It is noteworthy that even though the electronic temperatures reached 5.0 eV, the phonon frequencies of the IVB group carbides were positive. Thus the IVB group carbides have a higher stability under ultrafast laser irradiation, compared to metal (W, Na), semimetal (Bi) and semiconductor (Si, SiC).[9,35] All three of these become unstable in the electronic temperature range from 1.2 eV to 3.39 eV. The band gaps of the phonon density of states are 4.5801 THz, 2.624 THz, and 4.4339 THz for TiC, ZrC, and HfC at the ground state, respectively. Nevertheless, the phonon bands disappear at Te = 5.4 eV for TiC and Te = 5.17 eV for ZrC, and the phonon bands decrease to 1.3496 THz at Te = 5.5 eV for HfC.

Fig. 2. (color online) The phonon spectra of (a) TiC, (b) ZrC, and (c) HfC obtained as a function of electronic temperature Te.
3.3. Thermodynamic properties at different electronic temperatures

Based on the phonon frequencies and phonon density of states, we investigated the thermodynamic properties (the phonon free energy F in Fig. 3, the phonon internal energy E in Fig. 4, the phonon entropy S in Fig. 5, and the phonon heat capacity CV in Fig. 6) of ZrC at different electronic temperatures. These results provide reliable data for further theoretical and experimental studies on the relevant properties of the transition metal carbides. In Fig. 3, we can see that the phonon free energies increase as the lattice temperature is elevated. Furthermore, the phonon free energy–lattice temperature curves shift uniformly down with an increase in the electronic temperature. From Te = 0 eV to 1 eV the relationship between the phonon free energy and the electronic temperature is almost unchanged. The FT curves move down gradually with the elevation of the electronic temperature. Figure 4 indicates that the electronic excitation negatively affects the phonon internal energy of ZrC when the lattice temperature is relatively low. However, the discrepancies of the phonon internal energy among the different electronic temperatures are not obvious, as the lattice temperature increases.

Fig. 3. (color online) The phonon free energy F of ZrC at different electronic temperatures.
Fig. 4. (color online) The phonon internal energy E of ZrC at different electronic temperatures.
Fig. 5. (color online) The entropy S of ZrC at different electronic temperatures.
Fig. 6. (color online) The phonon heat capacity CV of ZrC at different electronic temperatures.

The phonon entropies corresponding to the lattice temperatures are provided in Fig. 5 for various electronic temperatures. It can be seen that the phonon entropy S increases with the elevation of the lattice temperature at different electronic temperatures. The electronic excitation positively affects the phonon entropy S at the same lattice temperature, which confirms that the randomness of ZrC increases as the effect of electronic excitation is enhanced. It is noteworthy that the phonon entropy-temperature curves exhibit a jump-change between Te = 4 eV to 5 eV. Feng et al. proposed that such an abrupt change may be an indication of the target materials undergoing a phase transition, because of the imaginary transverse acoustic phonon frequencies in the jump-change.[45] However, there are no imaginary frequencies between Te = 4 eV to 5 eV in Fig. 2. Therefore, the viewpoint of Feng et al. is not suited to our calculations. We consider that the difference is associated with the characteristic bonding forms of ZrC.

The phonon heat capacity CV curves are provided in Fig. 6 at different electronic temperatures. At the low lattice temperature (less than 400 K), the effects of ultrafast laser irradiation on the phonon heat capacity are distinct. When the lattice temperature is larger than 900 K, the influence of ultrafast laser irradiation can essentially be neglected.

3.4. Elastic properties and the melting temperatures at different electronic temperatures

The elastic properties of a lattice are directly described by the elastic constants Cij. For the cubic lattice, there are three independent elastic constants (C11, C12, and C44).

It is possible to show that a cubic lattice is dynamically stable if , , and , criteria called the Born stability criteria. We defined the bulk modulus as and the elastic shear constant .[41] The computed quantities of elastic modulus and melting temperature Tm of TiC in Figs. 7(a) and 7(b), ZrC in Figs. 7(c) and 7(d) and HfC in Figs. 7(e) and 7(f) at different electronic temperatures are exhibited. We can see that both TiC and ZrC have approximately the same decreasing trends, except that the values are almost unchanged from Te = 0 eV to 1 eV. HfC obviously moves down until the electronic temperature reaches 4 eV. The elastic shear constants of the IVB group carbides satisfy the Born stability criteria at Te = 5 eV, which indicate the stability of the crystal for IVB group carbides.

Fig. 7. (color online) Electronic temperature-dependent elastic moduli, Debye temperature ΘD, and melting temperature Tm in (a) and (b) TiC, (c) and (d) ZrC, and (e) and (f) HfC.

When the target materials are in the ground state, the Debye temperatures of 920 K for TiC, 705 K for ZrC, and 547 K for HfC are close to the previously reported theoretical values of 935 K, 699 K, and 544 K.[46] The predicted melting temperatures of TiC decreased from 3340 K at the ground state to 1159 K at Te = 5 eV. In addition, the predicted melting temperatures of ZrC exhibited a similar trend from 3698 K at the ground state to 1415 K at Te = 5 eV. Finally, the predicted melting temperatures of HfC decreased from 4201 K to 2241 K, but these values are almost unchanged from Te = 0 eV to 4 eV. The Debye temperatures of the IVB group carbides decreased at the ground state as the atomic number of the IVB group increased, but the melting temperatures reveal an inverse trend. The predicted melting temperatures of the IVB group carbides decreased as the electronic temperature increased. Particularly, the predicted melting temperatures of HfC are higher than TiC and ZrC at different electronic temperatures. Therefore, the use of HfC is advantageous in extreme high temperature environments.

4. Conclusion

In the present work, the effect of ultrafast laser irradiation on the lattice stability of ceramic materials (IVB group carbides) was investigated using density functional theory (DFT). The lattice parameters at different electronic temperatures for the different structures were calculated and the values of the lattice parameter for the ground state are in excellent agreement with experimental values. Phonon frequencies were obtained using inter-atomic force constants in the real space within DFPT at different electronic temperatures. The calculated phonon frequencies of TiC and ZrC for the ground state show good agreement with previously reported calculations and experimental values. The phonon frequencies of the IVB group carbides are positive even though the electronic temperature peaked at 5.0 eV. Thus IVB group carbides are more stable under ultrafast laser irradiation compared to the metal (W, Na), semimetal (Bi) and semiconductor (Si, SiC), which is a beneficial property for nuclear and aeronautical applications. The thermodynamic properties of ZrC were calculated as a function of temperature for different electronic temperatures. We demonstrated that the elastic shear constants of the IVB group carbides satisfy the Born stability criteria at Te = 5 eV, which indicate the stability of the crystal for IVB group carbides. The Debye temperatures of the IVB group carbides decreased at the ground state as the atomic number of the group increased, but the melting temperatures exhibited an inverse trend. In addition, a comparison of the predicted melting temperature of the IVB group carbides revealed that HfC is best suited for extreme high temperature environment applications.

Reference
[1] Harrison R W Lee W E 2016 Adv. Appl. Cream. 115 294
[2] Hwu H Chen J G 2005 Chem. Rev. 105 185
[3] Price D L Cooper B R 1989 Phys. Rev. 39 4945
[4] Craciun D Socol G Lambers E McCumiskey E J Taylor C R Martin C Argibay N Craciun V 2015 Appl. Surf. Sci. 352 28
[5] Upadhya K Yang J M Hoffman W P 1997 Am. Ceram. Soc. Bull. 76 51
[6] Reynolds G H Janivier J C Kaae J L Morlevat J P 1976 J. Nucl. Mater. 62 9
[7] Minato K Ogawa T Fukuda K Nabielek H Sekino H Nozawaw Y Takahashi L 1995 J. Nucl. Mater. 224 85
[8] Yutai K Gokul V Takashi N Lance L S 2013 J. Nucl. Mater. 441 718
[9] Shen Y H Gao T 2015 J. Alloy. Compd. 645 193
[10] Mayrhofer P H Mitterer C Hultman L Clemens H 2006 Progr. Mater. Sci. 51 1032
[11] Amir S G Omid G 2012 Aircr. Eng. Aerosp. Tec. 84 75
[12] Lance L S Yutai K Sosuke K 2010 J. Nucl. Mater. 399 200
[13] Hojou K Otsu H Furuno S Sasajima N lzui K 1996 J. Nucl. Mater. 239 279
[14] Gan J Yang Y Dickson C Allen T 2009 J. Nucl. Mater. 389 317
[15] Huang Y Maier B R Allen T R 2014 Nucl. Eng. Des. 277 55
[16] Hu X D Nori F 1997 Phys. Rev. Lett. 79 4605
[17] Kitagawa R Takebe H Morinaga K 2003 Appl. Phys. Lett. 82 3641
[18] Cavalleri A Tóth Cs Siders C W Squier J A Raksi F Forget P Kieffer J C 2001 Phys. Rev. Lett. 87 237401
[19] Zijlstra E S Huntemann N Garcia M E 2008 New J. Phys. 10 033010
[20] Diakhate M S Garcia M E 2009 Phys. Rev. 79 094117
[21] Jeschke H O Garcia M E Lenzner M Bonse J Krüger J Kautek W 2002 Appl. Surf. Sci. 197 839
[22] Rousse A Rischel C Fourmaux S Uschmann I Sebban S Grillon G Balcou P Forster E Geindre J P Audebert P Gauthier J C Hulin D 2001 Nature 410 65
[23] Johnson S L Heimann P A Lindenberg A M Jeschke H O Garcia M E Chang Z Lee R W Rehr J Falcone R W 2003 Phys. Rev. Lett. 91 157403
[24] Kresse G Joubert D 1999 Phys. Rev. B 59 1758
[25] Kresse G Furthmuller J 1996 Phys. Rev. 54 11169
[26] Kohn W Sham L J 1965 Phys. Rev. 140 A1133
[27] Blochl P E 1994 Phys. Rev. 50 17953
[28] Perdew J P Burke K Ernzerhof M 1996 Phys. Rev. Lett. 77 3865
[29] Perdew J P Burke K Wang Y 1996 Phys. Rev. 54 16533
[30] Kresse G Furthmuller J Hafner J 1995 Europhys. Lett. 32 729
[31] Parlinski K Li Z Q Kawazoe Y 1997 Phys. Rev. Lett. 78 4063
[32] Baroni S Gironcoli S de Corso A D Giannozzi P 2001 Rev. Mod. Phys. 73 515
[33] Lee C Gonze X 1995 Phys. Rev. 51 8610
[34] Feng S Q Zang H P Wang Y Q Cheng X L Yue J S 2016 Chin. Phys. 25 016701
[35] Yan G Q Cheng X L Zhang H Zhu Z Ren D H 2016 Phys. Rev. 93 214302
[36] Souadkia M Bennecer B Kalarasse F 2013 J. Phys. Chem. Solids 74 1615
[37] Recoules V Clérouin J Zérah G Anglade P M Mazevet S 2006 Phys. Rev. Lett. 96 055503
[38] Nartowski A M Parkin I P MacKenzie M Craven A J MacLeod I 1999 J. Mater. Chem. 9 1275
[39] Yvelin G Szymon L D Dorothy M D Alexander L S Katsumi T 2014 Phys. Rev. 90 094103
[40] Mazevet S Clerouin J Recoules V Anglade P M Zerah G 2005 Phys. Rev. Lett. 95 085002
[41] Giret Y Gelle A Arnaud B 2011 Phys. Rev. Lett. 106 155503
[42] Pintschovius L Reichardt W Scheerer B 1978 J. Phys. C: Solid State Phys. 11 1557
[43] Weber W 1973 Phys. Rev. B 8 5083
[44] Isaev E I Ahuja R Simak S I Lichtenstein A I Vekilov Yu Kh Johansson B Abrikosov I A 2005 Phys. Rev. 72 064515
[45] Feng S Q Zhao J L Cheng X L Zhang H 2013 J. Appl. Phys. 114 043519
[46] Chauhan M Gupta D C 2013 Diam. Relat. Mater. 40 96